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We study multifield contributions to the scalar power spectrum in an ensemble of six-field infla- 
tionary models obtained in string theory. We identify examples in which inflation occurs by chance, 
near an approximate inflection point, and we compute the primordial perturbations numerically, both 
exactly and using an array of truncated models. The scalar mass spectrum and the number of fluctu- 
ating fields are accurately described by a simple random matrix model. During the approach to the 
inflection point, bending trajectories and violations of slow roll are commonplace, and 'many-field' 
effects, in which three or more fields influence the perturbations, are often important. However, in a 
large fraction of models consistent with constraints on the tilt the signatures of multifield evolution 
occur on unobservably large scales. Our scenario is a concrete microphysical realization of quasi-single- 
field inflation, with scalar masses of order H, but the cubic and quartic couplings are typically too 
small to produce detectable non-Gaussianity. We argue that our results are characteristic of a broader 
class of models arising from multifield potentials that are natural in the Wilsonian sense. 
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1 Introduction 



Inflation [U El E] provides a superb explanation for the observed spectrum of cosmic microwave 
background (CMB) anisotropies. The simplest and best-understood models of inflation involve 
a single field slowly rolling down a relatively flat potential, but more complicated models involv- 
ing multiple light fields and/or violations of slow roll are arguably more natural, and provide 
the prospect of distinctive signatures. 

Strong theoretical arguments motivate the consideration of inflationary models with many 
light fields. In theories with spontaneously broken supersymmetry, naturalness suggests that 
scalars receive masses that are at least of order the gravitino mass. If the inflationary energy 
is the dominant source of supersymmetry breaking in the early universe, scalars that couple 
with gravitational strength to this energy will acquire masses of order the inflationary Hubble 
parameter, H. Thus, moduli in the theory do not decouple from the inflationary dynamics, 
and can be light enough to fluctuate]^] This picture is well-attested in flux compactifications of 
string theory, which typically include tens or hundreds of moduli with masses clustered around 
H. In cases where the moduli potential is computable, one generally finds a complicated, high- 
dimensional potential energy landscape with structure dictated by the spectrum of Planck- 
suppressed operators in the theory. The nature of inflation in such a potential is an important 
and urgent problem. 

However, despite intensive efforts to understand multifield inflation over the past decade, 
most analyses of explicit models consider only two light fields. Moreover, even though significant 
analytical tools have been developed to trace the evolution of primordial perturbations outside 



the Hubble radius in models violating the slow roll approximation (see £2.3.2 for a brief review 
of prior results), a large majority of works on the subject do make a slow roll expansion during 
Hubble exit. An understanding of truly general models involving several fields with masses of 
order H and arbitrary dynamics remains necessary. 

In this work we study the primordial perturbations produced by inflation in a class of six- 
field potentials obtained in string theory [5], corresponding to a D3-brane moving in a conifold 
region of a stabilized compactificationj^] Our approach is statistical: instead of directly fine- 
tuning the potential to achieve inflation, we draw potentials at random from a well-specified 
ensemble and study realizations that inflate by chance. We then compute the exact dynamics of 
the linear perturbations numerically, making no slow roll approximation. To reveal the physical 
processes underlying the resulting perturbations, we recompute the perturbations using a range 
of approximate, truncated descriptions that retain different subsets of the entropic modes, and 
we then cross-correlate the exact and approximate answers. For example, we identify multiple- 
field contributions to the perturbations by comparing the exact spectrum to the spectrum 
calculated in a single-field truncation (in which no slow roll approximation is made). 

The organization of this paper is as follows. In §|2] we describe our method for generating 



^ee [4] for a discussion of some of the cosmological signatures of models with spontaneously broken super- 
symmetry. 



2 This system was recently studied by Dias, Frazer, and Liddle in [6]. In ^ 4.1 we explain how our conclusions 
are qualitatively different from those of [5]: in contrast to [BJ, we find that an adiabatic limit is reached during 
inflation, so that the curvature perturbations can be predicted without modeling the details of reheating. 
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inflationary trajectories and computing the corresponding primordial perturbations. In £[3] we 
show that the scalar mass spectrum follows from a simple matrix model, and we assess the 
incidence and consequences of slow roll violations. In §|4] we study two-field and many-field 
contributions to the scalar power spectrum. We conclude in §|5j 



2 Method 

In this section we briefly review warped D-brane inflation, describe how we identify an ensemble 
of inflating solutions, and then explain how we compute the primordial perturbations. 



2.1 Background evolution in D-brane inflation 

In warped D-brane inflation [7], inflation is driven by a D3-brane moving toward an anti-D3- 
brane in a warped throat region of a flux compactification. The structure of the potential 
for this configuration has been derived in [HI |5]. The operators in the effective Lagrangian are 
dictated by the throat geometry, and have been computed explicitly, but the Wilson coefficients 
are determined by the detailed configuration of sources (fluxes, Euclidean D-branes, etc.) in 
the bulk of the compactification. Thus, with present knowledge the Wilson coefficients can at 
best be modeled statistically. 

In [9], two of us, in collaboration with N. Agarwal and R. Bean, studied a large number 
of realizations of the potential, drawing the Wilson coefficients from statistical distributions, 
and found that although the typical scale of all six scalar masses is H, accidental cancellations 
among many terms could nevertheless lead to prolonged inflation. The inflationary trajectory 
took a characteristic form: the D3-brane initially moved rapidly in the angular directions of 
the conifold, spiraled down to an inflection point in the potential, and then settled into an 
inflating phase. It was established in jU] that the inflationary phenomenology has negligible 
dependence on the detailed form of the statistical distribution of the Wilson coefficients: a sort 
of universality emerges in this complicated ensemble]^] 

Our method for studying the evolution of the homogeneous background is identical to that 
of [9], to which we refer for further details. We began with an ensemble of more than 13 million 
potentials describing a D3-brane in a conifold geometry. For each potential we started with 
zero kinetic energy at a fixed location in the conifolcQ and evolved the background equations 
of motion numerically, as detailed in [9j. After discarding trials in which the inflaton became 
stuck in a local minimum or was ejected from the throat region, we obtained 18731 realizations 

3 See[ini[IIJ[Il[IS[Tl[IIl[M[^ related work on multifield evolution in D-brane 

inflation. 

4 Trajectories with nonvanishing initial kinetic energy were studied in [S], where it was found that as long 
as the initial kinetic energy is no larger than the initial potential energy, the effect on the phenomenology is 
minimal. Similarly, it was shown in [9] that the dependence on the initial radial location is minimal, while the 
rotational invariance of our ensemble of potentials implies that we do not lose any generality by fixing the initial 
angular positions. 
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yielding at least 66 e-fold^]of expansion followed by a hybrid exitj^] 

A word of caution about choices of measure is necessary. As shown in [9], the inflation- 
ary phenomenology is substantially independent of the statistical distribution from which the 
Wilson coefficients are drawn, and of the measure taken on the space of homogeneous initial 
conditions. However, we do not include a measure factor that weights histories according to 
the amount of expansion, and this must be borne in mind when interpreting our results. 



2.2 Aspects of inflection point inflation 

When inflation arises in the ensemble of potentials considered in this work, it does so at an 
approximate inflection point [23l EU |25j |9] . Before proceeding we will recall a few elementary 
properties of single-field inflection point inflation that have significant ramifications for our 
analysis. 

Inflection point inflation begins in a region of field space where the potential is positively 
curved and evolves to a region where the curvature is negative (see Fig. [I}. As the size of the 
inflection point region in Planck units is typically extremely small (see e.g. [9]), the potential 
slow roll parameters ey = (^) and t]v = M^^y- obey ey <C \vv\, so that n s ~ 1 + 
2r]v- Hence, the scalar power spectrum calculated in the single-field slow roll approximation 
is initially blue and becomes red. Correspondingly, the tilt measured in the CMB is dictated 
by the point in field space at which observable modes exit the Hubble radius. Inflection points 
producing a total of N e f» 120 e-folds of inflation have approximately 60 e-folds below the 
inflection point and 60 e-folds above, so that the primordial spectrum on large angular scales is 
scale-invariant. Inflection points producing N e > 120 e-folds have red spectra, while inflection 
points producing N e < 120 e-folds have blue spectra and are in tension with observations. 



2.3 Perturbations in multifield inflation 

In this section, we review an efficient approach to computing the primordial perturbations in 
multifield inflation, following [26] and |27j . 



2.3.1 Equations of motion for exact multifield treatment 

Our starting point is the action for a collection of scalar fields (ft 1 , endowed with a metric^ 
Gu(4> K ) on field space, interacting through a potential V(0 / ), and minimally coupled to gravity: 




5 Although we are interested in requiring at least 60 e-folds of expansion to solve the horizon problem, we 
insist on a total of 66 e-folds in order to be able to impose Bunch-Davies initial conditions well before the CMB 
exits the Hubble radius. 

6 In practice we take the end of inflation to occur when the D3-brane reaches a fixed radial location slightly 
above the tip. 

7 Relativistic motion, corresponding to Dirac-Born-Infeld inflation [2§], did not arise in our ensemble 
[501 0|] , so it suffices to consider the two-derivative kinetic term. 
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Figure 1: A single-field inflection point potential, taken from [25]. If the inflaton is above the 
inflection point 60 e-folds before the end of inflation, V" > and the scalar power spectrum 
is blue, as indicated by the shading. A red spectrum requires that the inflaton has passed the 
inflection point 60 or more e-folds before the end of inflation. 



(see [31] for the explicit expression for the field space metric in the case of the conifold). 

The background metric is assumed to be of the spatially flat Friedmann-Lemaitre- Robertson- 
Walker form 

ds 2 = -dt 2 + a 2 {t)dx 2 , (2) 

where t is cosmic time and a(t) denotes the scale factor. If the fields depend only on t, the 
equations of motion take the simple form 

3H 2 = \<r 2 + V, (3) 

H = Af\ (4) 
Vtft + 3H^ + G IJ Vj = 0, (5) 

where dots denote derivatives with respect to t, H = a/a is the Hubble parameter, \a 2 = 
^Gij(f> I <p J is the kinetic energy of the fields, and, here and in the following, T>tA l = A 1 + 
Tj K (f) J A K for a field space vector A 1 . 

The dynamics of linear perturbations about this background is governed by the second-order 
action [321 EQ EE 

S {2) = j dtd 3 xa 3 \G IJ V t Q I V t Q J - ^Gjjd^d'Q- 7 - M 7J Q 7 Q J ) , (6) 
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where the Q 1 are the field fluctuations in the spatially flat gauge and the mass (squared) matrix 
is given by 

(7) 



M u = Vjj - K IKLJ j> K 4> L - \v t 



Here Vjj = Vjj — YfjV t K is the covariant Hessian, IZiklj is the Riemann tensor associated to 
the field space metric, and field space indices are raised and lowered using Gjj. From equation 
^ one easily deduces the equations of motion for the linear fluctuations (in Fourier space): 

VtVtQ 1 + mVtQ 1 + ^Q 1 + M]Q J = . (8) 

(X 

2.3.2 The adiabatic/entropic decomposition 

Following [33j EE], it is useful to decompose the field fluctuations into the so-called (instanta- 
neous) adiabatic and entropic perturbations. The adiabatic perturbation is defined as 

Q a = e^Q 1 , (9) 

where e£ = <fi ! /& is the unit vector pointing along the background trajectory in field space, 
whereas entropic fluctuations represent fluctuations off the background trajectory. The adia- 
batic fluctuation is directly proportional to the comoving curvature perturbation 71, 

K = —Qa , (10) 

a 

while the genuinely multifield effects are embodied by the entropic fluctuations. 

One of the entropic modes plays a distinguished role: this is the 'first' entropic fluctuation 
Q S1 = e sl iQ I , which is the fluctuation along the direction of acceleration perpendicular to the 
background trajectory, where 

ej=- /"^ =, (11) 
a/ -L /j VjVj 

and _L /J = G IJ — e^e^ is the projection operator on the entropic subspace. The first entropic 
mode instantaneously couples to the adiabatic perturbation, but the remaining entropic modes 
do not. 

The adiabatic equation of motion can be written in the compact form 

Q a + 3HQ a + ^ + m 2 )j Q a = {2Ht 1i _Q Si \ -f^ + YA 2H V± Q Sl , (12) 

where 

V, - (13) 

is a very important dimensionless parameter measuring the size of the coupling between the 
adiabatic mode and the first entropic mode. Here V Sl = e^Vj, and similarly for analogous 
quantities, and the adiabatic mass (squared) m 2 a is given by 

9 

TfL 

-^ = ^- rj 2 + t7 1 + 7 1 /H, (14) 
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with e = — -fjz an d V = ~\w- 

A brief overview of existing methods to study cosmological fluctuations in multifield infla- 
tion is appropriate at this stage. Several analytical methods have been developed to follow the 
evolution of perturbations outside the Hubble radius without making any slow roll approxi- 
mation, including the SN formalism [Ml |32| |35| [36] , the transfer functions method [37], the 
gradient expansion method [3E1 EH1 SB HI] , the covariant method [1211131111], and the transport 
equations method (321 HHJ ETJ HE] . A number of works incorporate perturbative corrections to 
the slow roll approximation during the epoch of Hubble exit, including [261 SHI EH ED E21 [53] . 
but exact results are scarce for systems with more than two fluctuating fields. Related investi- 
gations of the effects of heavy fields in multifield inflation include [5U [551 ESI E3 E3 EH EQl [61] , 
while numerical studies for two-field systems include [62] EH E3] ■ 

We stress that an exact numerical treatment is needed, at present, for models involving 
significant bending around Hubble crossing, which as we will see are common in our ensemble. 

2.3.3 Methods for computing the perturbations 

For each of the 18731 realizations of inflation yielding at least 66 e-folds of expansion followed 
by a hybrid exit, we numerically integrated the exact equations of motion ^ for the linearized 
perturbations corresponding to the scale exiting the Hubble radius 60 e-folds before the end 
of inflation. To relate this result to conceptually simpler models, we also computed the per- 
turbations using an array of approximate descriptions, which we now specify. To distinguish 
references to these specific models from more general uses of words such as "exact", "naive", 
etc., we will put the model name in a distinctive font. 

The naive model simply computes at Hubble exitj^] and incorporates neither slow roll 
violations nor multifield effects. The adiabatic or one-field model makes no slow roll approxima- 
tion, but discards the effects of all entropic modes: only fluctuations tangent to the trajectory 
at any given time are retained. 

The two-field model keeps the instantaneous adiabatic and first entropic fluctuations. For 
this case the first entropic equation of motion takes the form 

Q n + 3HQ S1 +(^ + m 2 s )j Q S1 = -2& V± n , (15) 

where 

m 2 si = V ;S1S1 - a 2 R Sluasi - H 2 rj]_ . (16) 

The three-field, four-field, and five-field models (which we will not use in this work) retain 
additional entropic perturbations along corresponding basis vectors in the decomposition of 
[26J. We will refer to the six-field model, which includes all five entropic modes and therefore 
incorporates all the physics of the linear perturbations, as the exact model. The scalar power 
spectra resulting in each model are named in a similar manner, e.g. V orie is the scalar power 
computed in the one-field model. 

8 All quantities described as 'evaluated at Hubble exit' in this work are actually averaged from one e-fold 
before Hubble exit until one e-fold after Hubble exit, in order to smooth quantities that may be rapidly changing. 
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By comparing the results of these approximate descriptions, one can unambiguously identify 
certain physical effects in the perturbations. Specifically, V one ^ "Pnaive signifies violations of slow 
roll, while V exact 7^ Vu, demonstrates that at least k + 1 fields contributed to the perturbations. 

Finally, a discussion about numerical implementation of the quantization and evolution of 
the perturbations is in order. As is well known (see for instance [64, 65, 62J), to determine 
the late-time power spectrum by numerically evolving the coupled equations ([8]), one cannot 
simply impose a single choice of initial conditions and solve the system ^ once to obtain the 
six perturbations Q 1 . This would introduce spurious interference terms in the power spectrum 
between what are supposed to be independent variables. On the contrary, one should identify 
six variables that are independent deep inside the Hubble radius, each corresponding to an 
independent set of creation and annihilation operators, and solve the system of equations ^ six 
timesj^jeach time imposing the Bunch-Davies initial conditions for only one of the independent 
variables, while setting the other variables to zero initially]^] One then extracts power spectra 
by summing the relevant quantities over all six runs. This is the numerical analogue of the fact 
that the various creation and annihilation operators are independent, so that their effects add 
incoherently (see £4.3.4 for a precise illustration of this method). 

Deep inside the Hubble radius, one can neglect the mass matrix in equation (|6]), so that 
identifying a set of independent variables is equivalent to identifying a set of vielbeins for the 
field space metric Gjj, which is easily accomplished numerically. We follow this strategy for the 
quantization, imposing initial conditions six e-folds before Hubble crossing, while still solving 
the system of equations ^ in the natural coordinate basis on the conifold (i.e. r, 6\, (pi, 9 2 , 02, ip 
— see [31]), which we find numerically efficient. 

The method above applies to the exact model, but some modifications are required for 
the k-field models with k < 6, where by definition only k independent variables are included. 
While there is no possible subtlety for the one-field model, for the two-field model for instance, 
the easiest set of vielbeins consists of the adiabatic and first entropic basis vectors. From a 
numerical perspective, the kinematical basis of [26] corresponds to a particular set of vielbeins 
with transparent physical meaning. 



2.3.4 Terminology for restricted ensembles 

For convenience, we now define important subsets of our ensemble. The full ensemble consists 
of all realizations yielding at least 66 e-folds of inflation. The effectively single-field subset 
consists of all realizations in which multifield contributions to the scalar power are at most 1% 
corrections: specifically, we require that 

£exact/one = \V exact /V one - l| < 0.01 . (17) 

9 Equivalently, one can solve a matrix- valued system once — see |5SJ EH EI]- We thank Ian Huston for 
pointing out Ref. |66) . 

10 In the k-field models with k < 6 described above, we set 6 — k entropic modes to zero throughout the 
evolution, not just initially. 
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The effectively multifield ensemble is the complementary subset, consisting of all cases with 
^exact/one > 0.01. Similarly, we define 

exact/two = Inexact /^two 1 I > (^"^0 

which measures the extent to which a two-field description is insufficient. The set of models 
with ^exact/two > 0.01 may be termed 'effectively many-field'. Finally, the observationally allowed 
ensemble consists of all realizations satisfying the WMAP7 constraints on the tiltp]of the scalar 
power spectrum at 2a. 



3 The scalar mass spectrum and violations of slow roll 

Before describing the primordial perturbations, we will first characterize the spectrum of scalar 
masses. We begin with empirical observations about the mass spectrum (§|3.1 ), and then turn to 



obtaining several of these properties from a random matrix model ({ 3.2 ). In £ 3.3 we present key 
consequences of the scalar mass spectrum, focusing on violations of slow roll and the resulting 
imprint in the power spectrum. 

3.1 Properties of the mass spectrum 

For each inflationary trajectory we computed the Hessian matrix V'j of the scalar potential, 
in terms of the canonical coordinates constructed at the relevant location on the conifold, at 
the moment that the CMB exited the Hubble radius. We denote the ordered eigenvalues of 



the Hessian as m{ < . . . < m^, with corresponding eigenvectors ip^ . . . ip 6 . A histogram of 



m\ < . . . < m\ appears in Fig. [2] 

The mass spectrum of ipi has a sizable spike at mildly tachyonic values: m\ > — O.lif 2 . The 
lower bound on the mass-squared is a consequence of conditioning on prolonged inflation: at 
a generic point in the field space, not along an inflationary trajectory, ip\ would have a much 
stronger tachyonic instability. In contrast, ip2 is almost never tachyonic: we found m\ < in 
only 4 out of 18731 realizations. Moreover, ip 3 through ip Q were never tachyonic in our ensemble. 

An important property of the spectrum is that m\ ~ m\ and m\ ~ wig, to an accuracy 
of a few percent. This is not a statistical statement: these eigenvalues are degenerate in each 
realization, not just after taking the ensemble average. The corresponding histograms are 
overlaid in Fig. [2] 

Although we have evaluated the masses at Hubble exit, each mass slowly changes during 
the course of inflation, in a predictable way. We find that the four heaviest fields become 
fractionally more massive: in one e-fold around Hubble crossing, there is a change <5m 2 ~ .02m 2 
for a = 3 ... 6. In contrast, the two lightest fields diminish in mass, with 6m\ ~ — .01m 2 . and 
5m\ ~ —0.181m!! 2 . 



for a discussion of constraints on the amplitude and running of the scalar power spectrum. 



n See [4.4 

12 In general the ipi are complicated combinations of the six natural coordinates on the conifold, though in 
trajectories that are primarily radial there is significant overlap between ipi and the radial coordinate r. 
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Figure 2: The mass spectra of the six scalar fields, in units of H 2 . The leftmost peak, which 
has support at tachyonic values, corresponds to the lightest (adiabatic) field ipi- The next peak 
corresponds to the second-lightest field ip2- The third peak corresponds to ips and ip4, which 
are nearly degenerate in each realization, and the broad final peak similarly corresponds to ips 
and ip 6 . 

The masses that are relevant for determining the evolution of the fluctuations consist of 
a standard part coming from the potential, namely the V-u term in equation ([7]), as well 
as kinematic contributions involving time derivatives of background quantities, corresponding 
to the remaining terms in equation ([7]). We find that the kinematic contributions introduce 
corrections at the level of one part in 10 4 , and so can be neglected for practical purposes. 

3.2 A random matrix model for the masses 

We will now show that distinctive qualitative features of the mass spectrum can be understood 
using random matrix theory. The framework for this analysis is 'random supergravity', by 
which we mean an ensemble of four-dimensional Af = 1 supergravity theories whose Kahler 
potential and superpotential are random functions, in a sense made precise in [68J. A matrix 
model that is slightly simpler than that of [68J will suffice for our purposes: we take the Hessian 
matrix % to be of the formF*! 



13 Cf. [68 for an explanation of the relationship between this model and the full Hessian matrix in supergravity. 
We thank T. Wrase for helpful discussions of this point. 




(19) 
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up to a shift proportional to the identity matrix, where A, B, and C are 3x3 complex symmetric 
matrices. We take the entries of A, B, and C to be random complex numbers drawn from a 
normal distribution; as explained at length in [68J , the spectrum of ~H is essentially independent 
of the statistical properties of the matrix entries. For C = the spectrum is positive-definite 
and doubly degenerate, but C ^ breaks the degeneracies and permits negative eigenvalues. 

By requiring inflation, one has effectively imposed a lower bound on m\ , which clearly affects 
the empirical spectrum shown in Fig. |2j We should therefore impose a corresponding restriction 
in the matrix model: we take m\ > —OAH 2 . (The overall scale is arbitrary in the matrix model, 
and we set our units by matching the right tail of the rightmost peak in Fig. |2j) In Fig. [3] we 
show the results of simulations of the spectrum in this simple matrix model. By adjusting 
parameters — such as the weights assigned to A, B, and C, or the relative variance of their 
entries — one can achieve reasonably good quantitative modeling of the empirical spectrum, 
but it is not clear that this is physically meaningful. We have instead presented the results 
for the simplest case, with equal weights for all matrices, and equal variance for all entries, 
for which the qualitative agreement is already surprisingly good. Notice that the eigenvalues 
are approximately pairwise degenerate, which matches the empirical result for m\ w m\ and 
m\ ~ wig but is a less accurate model of m\ and mf. In contrast to the analysis of [68], the 
matrices in question are not large: they are 3 x 3, so that large N arguments (where N is the 
size of the matrix) are marginal at best. However, for C = and any N, the spectrum of H is 
exactly doubly degenerate, and the approximate degeneracies in the mass matrix with C ^ 
hold at finite N. 
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Figure 3: The mass spectrum simulated in the matrix model given in equation (19), cf. [68J. 
Comparing to Fig. [2j we see that the model is in good qualitative agreement with the results 
of simulations in the ensemble of inflaton potentials. 



Another key feature of the empirical spectrum that is reproduced in the matrix model is 
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the sharp left edge in the probability density for m\. The requirement of N e > 66 e-folds of 
inflation effectively imposes a lower bound m\ > — 0.1H 2 . The steep pileup of the probability 
density near this edge is a characteristic feature of eigenvalue spectra of ensembles of random 
matrices subject to constraints on the smallest eigenvalues. 

We conclude that qualitative features of the mass spectrum arising from the scalar potential 
on the conifold can be reproduced in the very general random matrix model of [68], or its simpli- 



fied version (19), which governs any M = 1 supergravity theory for which the Kahler potential 
and superpotential are random functions to good approximation. Correspondingly, these fea- 
tures are plausibly properties of a general random supergravity theory, not consequences of the 
particular conifold setting of the present work. This strongly motivates using kindred matrix 
models to study much more general many-field inflationary scenarios. 

To make a further observation about the scope of the matrix model, we briefly recall the 
method used in [3] to compute the scalar potential for a D3-brane on the conifold. The com- 
putation of [5] amounted to determining the most general solutions of ten-dimensional super- 
gravity, in expansion around the Klebanov-Strassler solution, with certain asymptotics. As 
such, these solutions incorporate in full detail the structure of a conifold geometry attached to 
a stabilized compactification, but — being intrinsically ten-dimensional — lack any manifest 
four-dimensional M = 1 supersymmetry. In particular, the Kahler potential and superpotential 
of the four-dimensional theory arising upon dimensional reduction were not obtained directly 
in [5]. We therefore find it remarkable that a random matrix model based on four-dimensional 
random supergravity is in excellent agreement with our calculation of the mass spectrum arising 
from the ten-dimensional results of [5]. 

A further quantity of significant interest is the number of scalar fields that are light enough 
to fluctuate during inflation. We will now show that an efficient measure of the number of 
fluctuating scalars is a weighted average that takes into account the reduced contributions of 
fields with masses m — > 3/2H that barely fluctuate. To determine the proper weighting, we 
consider the equation of motion for the perturbations. The canonically normalized perturbation 
v = aQ of a test scalar field of mass m in de Sitter spacetime obeys the equation of motion, in 
Fourier space, 

4+(k 2 -±-Jv 2 -]))v k = 0, (20) 



where 



o 9 m 2 , . 

" S 4"^' < 21) 

and ' denotes a derivative with respect to conformal time r. The solution of this equation with 
the appropriate Bunch-Davies behavior inside the Hubble radius reads 

v k (r) = ^e^/^^H^i-kr) , (22) 

where is the Hankel function of the first kind, of order v, and we use the convention 

v = when z/ 2 > , (23) 

v = %\/-v 2 when v 2 < . (24) 
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The power spectrum of fluctuations at Hubble crossing, i.e. for — kr± = 1, can then be computed 

as 



V, 



Q k ,m 



2 T* 



2tt ; 



:\Qk(r* 



(25) 



For each of our inflationary realizations, one can sum the six contributions of the form (25) 



corresponding to the six eigenvalues of the Hessian matrix determined in {3A_ We thus obtain 
an analytical estimate of the total power spectrum of scalar field fluctuations at Hubble exit, 



(26) 



In realizations with small bending of the trajectory (see £4.2 ), the estimate (26 ) agrees extremely 



well with a direct numerical calculation of the total power spectrum of field fluctuations at 



Hubble exit, with a precision of order 0.01%. In more general models, the estimate (26) is only 



qualitatively correct, and our full numerical treatment is necessary (see for instance ^4.3.4 ). 
The effective number of fluctuating fields, n/, can then be defined by comparing the estimate 



(26) to the corresponding expression for a massless scalar field: 



(27) 



In Fig. [4] we show a histogram of n/ in our ensemble of inflationary realizations. Remarkably, 
the histogram is sharply peaked, and nf falls between 2 and 3 in about 99% of our realizations. 

We caution the reader that the number of fields that are light enough to fluctuate, nj, is in 
general different from the number of fields n^. that contribute to the curvature perturbation: 
for example, if the background trajectory is straight, then entropic perturbations are not con- 
verted to curvature perturbations, and n-ji = 1 for any nj. The calculation above concerns the 
mass spectrum, and hence determines the statistical distribution of n/. We will determine the 



distribution of n-ji in £4.3.3 



3.3 Violations of slow roll 

A pivotal property of the scalar mass spectrum discussed in the preceding section is that the 
lightest field ipi has a mass m\ ~ H 2 in a large fraction of realizations yielding N e > 66 e-folds 
of inflation. Thus, the slow roll approximation is only marginally applicable. 

In Fig. [5] we show that the mass-squared of the adiabatic direction, evaluated 60 e-folds be- 
fore the end of inflation, has a clear dependence on the total number of e-folds, N e . Realizations 



with N e < 100 have ml/ H 2 > 1, while realizations with N e > 100 havej^j^a ~ -0.1H 2 . This 
is simply a consequence of the fact that inflation is occurring at an inflection point. Realiza- 
tions yielding N e < 120 have the 60 e-fold mark above the inflection point, where the curvature 
of the potential is positive, while realizations yielding more inflation have the 60 e-fold mark 



slightly below the inflection point. See £2.2 



14 It is straightforward to show analytically that 60 e-folds before the end of inflation in a single-field inflection 
point model yielding N e 3> fOO, the inflaton mass obeys m 2 w — in excellent agreement with our simulations. 
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Figure 4: Histogram showing the relative probability of the number n/ of scalar fields that are 



light enough to fluctuate during inflation, cf. equation (27) 




Figure 5: The mass-squared of the adiabatic fluctuation in units of H, evaluated 60 e-folds 
before the end of inflation, versus the total number of e-folds, N e . 



Violations of slow roll provide significant corrections to the scalar power spectrum in our 
ensemble. In Fig. [6] we show the ratio of the one-field power to the naive power, as a function 
of m 2 /H 2 . We learn that increasing the mass of the adiabatic direction decreases the scalar 
power in a predictable manner. 

The effects of slow roll violations on the tilt for effectively single-field models are shown in 
Fig. [7j Except for the handful of cases with m 2 /H 2 > 1.3, the exact tilt is more blue than 
the naive tilt, by an amount that is strongly correlated with m 2 /H 2 . Thus, slow roll violations 
tend to shift the spectrum to be slightly more blue. 
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Figure 6: The ratio of the one-field power to the naive power, versus the adiabatic mass-squared 
in units of H. The left panel shows the entire ensemble, while the right panel is restricted to 



effectively single-field realizations (see £2.3.4). 
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Figure 7: The difference between the exact and naive spectral tilts, versus the adiabatic mass- 



squared in units of if, for effectively single-field realizations (see £2.3.4) 



4 Multifield effects 



We now turn to our primary objective, the characterization of multifield contributions to the 
scalar power spectrum. To begin, in §4.1 we discuss the degree to which multifield inflation in 
our ensemble is predictive. Next, in §4.2| we quantify the degree of bending of the trajectory. 



In £4.3 we characterize the frequency with which multifield and many-field effects arise, and in 



£4.4 we describe the consequences of imposing observational constraints on the tilt. We discuss 



the prospect of observable non-Gaussianities in £4.5 
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4.1 Decay of entropic perturbations 



Entropic perturbations are essential for super-Hubble evolution of the curvature perturbation: 
if all entropic modes become massive and decay at some point after Hubble exit, the curvature 
perturbation becomes constant, and one says that an adiabatic limit has been reached (see 
[69| W0\ I7T| 172] for recent discussions). Provided that an adiabatic limit is reached during 
the inflating phase, one can predict the curvature perturbation on observable angular scales 
without knowing the details of the end of inflation or of reheating. Conversely, failure to reach 
an adiabatic limit makes predictions contingent on an understanding of reheating. 

The approach to an adiabatic limit clearly depends on the masses of the entropic modes: 
modes with m > H decay quickly outside the Hubble radius, and those with m > \H do 
not oscillate at all. One might a priori expect that in our ensemble all six fields have masses 
of order H, with a substantial likelihood that more than one field is tachyonic. Our findings 



from £3.2 differ from this expectation in important details: there is at most 15 one tachyonic 
instability, the second-lightest field ip2 has m\ ~ H 2 , and the heavier fields ips ... ipe have masses 
considerably larger than H. Correspondingly, we expect that the entropic modes will decay 
during the course of inflation, so that an adiabatic limit is reached before the end of inflation. 

To see that this is the case, we examined a subset of realizations and verified that the total 
power of the five entropic perturbations decays exponentially. In 60 examples, the power in 
entropic modes at the end of inflation was never larger than 10~ 10 times the adiabatic power. 
In Fig. [8] we show the characteristic decaying behavior for 19 examples. 

One possible subtlety is worth mentioning. In the argument above we have assumed that 
the lightest field ipi corresponds to the adiabatic direction CT , so that having m\ ~ H 2 and 
m\, m\, ml, ml ^> m\ ensures that all five entropic modes are massive enough to decay quickly. 
There is good evidence for this identification: the histograms of m\ and of m 2 are nearly 
indistinguishable, but are quite different from those of the remaining fields. 

Somewhat different results concerning the approach to the adiabatic limit in a similar system 
have been reportecf^jby Dias, Frazer, and Liddle [B], who find that the angular fields have masses 
that are very small compared to H . Correspondingly, no adiabatic limit is reached, and there 
is a problematic persistence of entropic modes. Based on the analysis of [5] and [9], or more 
generally on the grounds of naturalness in a theory with spontaneously broken supersymmetry, 
we would expect masses of order H for all six fields (modulo a suppression of the mass of the 
adiabatic direction resulting from conditioning on prolonged inflation), which is consistent with 
the results reported here, but does not appear consistent with [5J. 



4.2 Bending of the trajectory 

Because the entropic perturbations generally decay rapidly outside the Hubble radius, they 
are relevant for the late-time curvature perturbation only if the inflationary trajectory bends 
shortly after they exit the Hubble radius, so that the entropic perturbations are rapidly con- 
verted to curvature perturbations. We therefore turn to characterizing the incidence of bending 

15 Strictly speaking, we found two tachyons in a negligible fraction of trials, in 4 out of 18731 examples. 
16 We thank Mafalda Dias and Jonathan Frazer for discussions of this point. 
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Figure 8: The curves show the rapid decay of the total power of the five entropic modes over 
the course of inflation, for 19 realizations. The decaying behavior depicted is representative 
of the ensemble: the lightest entropic mode generically has m 2 > H 2 , and was tachyonic at 
Hubble crossing in a negligible fraction of trials. 



trajectories in our ensemble. 

A very useful measure of turning is the parameter t]± defined in equation (13), which mea- 
sures the acceleration of the trajectory transverse to the instantaneous velocity. Analytic meth- 
ods that make a slow roll approximation around the time of Hubble crossing generally require 
i]± <C 1, which can be related to the 'slow-turn approximation' discussed in [S3]. Correspond- 
ingly, effects requiring r]± > 1 at Hubble crossing are only partially understood. Our treatment 
makes no approximation, and indeed r]± can be quite large: see Fig. |9j which shows the evo- 
lution of r]± in the first 15 e-folds of inflation in a handful of representative examples. The 
gradual decay of r]x can be understood from the properties of the evolution near the inflection 
point: while the inflaton is spiraling down to the inflection point, turning is generic, but in the 
immediate vicinity of the inflection point, where prolonged inflation occurs, the trajectory is 
quite straight. 

In the left panel of Fig. 10 we plot ?7_i_, evaluated at Hubble crossing, versus the total number 
of e-folds, N e . Evidently, T]± at Hubble crossing is quite small in scenarios giving rise to N e 3> 60 
e-folds of inflation, but is significant in scenarios yielding less inflation. This is consistent with 
the picture described above in which there are large transient contributions to rj± that decay 
after prolonged inflation. 

Another practical measure of the amount of turning is the 'total turn' AO, which we define 
to be the integral of rj± from six e-folds before Hubble crossing until the end of inflation: 
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Figure 9: The curves show the rapid changes in r]± during the first 15 e-folds for 9 realizations 
of inflation. Note that N e = refers to the start of inflation in each realization, rather than 
to the moment when the CMB crosses the Hubble radius, which occurs much later in many 
examples. 




Figure 10: Left panel: the slow roll parameter 7]± evaluated at Hubble crossing, versus the 
total number of e-folds, N e . Right panel: the total turn of the trajectory, from six e-folds 
before Hubble crossing until the end of inflation, A6 = J^_ G dN e T]±, versus the total number 



of e-folds, N P . 



Ad = lHC-edN eVx 
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From the right panel of Fig. 10 we see that the total turn is quite 



17 Turns occurring shortly before Hubble crossing can have an effect on the perturbations, and we therefore 
define the total turn as the integral of r]± from six e-folds before Hubble crossing to the end of inflation. 
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large in cases with N e < 80 e-folds of inflation. Notice that even in cases with arbitrarily 
many e-folds, a total turn of order AO ~ 1/2 remains likely. This is an effect from the end of 
inflation: as the inflaton falls off the inflection point, its trajectory very often bends (see [73] 
for a discussion of related points). We have verified this by checking that the restricted integral 
Jhc-6° dN e r)± is small in cases with N e ^> 60 e-folds, even though J^_ 6 dN e rj± is not. 



4.3 Multifield effects and many-field effects 

Having understood the properties of the mass matrix and of the background evolution that 
influence the evolution of entropic perturbations, we are in a position to understand multifield 
contributions to the primordial perturbations. 

One of the most interesting questions for a model of inflation with more than two fields 
is whether the primordial perturbations are effectively governed by only two fields, or instead 
have distinctive 'many-field' signatures (see [71] for a recent discussion). In a general iV-field 



model of inflation, as we reviewed in £2.3.2, one of the N — 1 entropic fluctuations, which 
we call the first entropic fluctuation, plays a distinguished role by instantaneously coupling to 
the adiabatic perturbation. Most explicit studies of inflation with more than one field have 
taken N = 2 for simplicity, in which case the entropic subspace is one-dimensional and this 
distinction is unnecessary. However, for our system there are five entropic modes, each of which 
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can in principle contribute to the curvature perturbation!^ We must therefore study not just 



the incidence of multifield effects, but also the incidence of 'many-field _ effects, in which the 
curvature perturbation is not dictated by the adiabatic and first entropic fluctuations alone, 
and instead receives contributions from the higher (second through fifth) entropic modes. The 



method for this analysis was described in £ 2.3.3 when V exact ^ V^, then at least k + 1 modes 
contribute to the curvature perturbation. 



4.3.1 Illustrative examples 



We begin by examining a collection of illustrative examples, shown in Fig. [TTJ For each example 
we display the power spectrum as a function of the number of e-folds since Hubble crossing, 
using the naive, one-field, two-field, and exact models. 

In the upper left panel we show a striking example for which all four models give different 
results. Around Hubble crossing the various models give roughly comparable results for the cur- 
vature perturbation, but conversion of entropic perturbations to curvature perturbations after 
Hubble crossing causes the two-field and exact results to grow sharply. As P eX act ^ V ne, multi- 
field effects cannot be ignored, but because also V exact ^ V two , we see that a two-field description 

18 Although the second through fifth entropic modes do not couple instantaneously to the adiabatic pertur- 
bation, they do couple to the first entropic mode, which can then source the adiabatic perturbation when the 
trajectory bends. 

19 For the purpose of this discussion, because 'multifield inflation' means 'inflation with two or more light 
fields', we will use 'many-field inflation' to refer to 'inflation with three or more light fields'. An important and 
challenging problem is to understand inflation with N 3> 1 light fields, but we do not use 'many' in this sense 
in the present work. 
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is inadequate. Notice the scale: many-field contributions to the curvature perturbation dwarf 
the two-field contribution, which is itself non-negligible. 

Next, in the upper right panel of Fig. 11 we exhibit a case in which V exact ~ T^two 7^ "^one, 
so that two-field effects cannot be omitted, but many-field effects are negligible. The lower 
left panel shows an example in which Pexact 7^ V^o ~ V one , for which many- field effects are 
significant but two-field effects are negligible. Finally, the lower right panel shows a rare example 
in which V exact ~ V two 7^ V one , so that two-field but not many-field effects are important, and 
the resulting power spectrum is compatible with observational constraints on the tilt. 





Figure 11: Scalar power spectra as functions of the number of e-folds since Hubble crossing, 
for four different realizations. The purple, orange, blue, and red lines correspond to the naive, 
one-field, two-field, and exact spectra, respectively. Lines with combined colors indicate that the 
corresponding spectra overlap to high accuracy. Upper left: many-field and two-field effects are 
both important. Upper right: two-field effects are present, but many-field effects are negligible. 
Lower right: two-field effects are present but many-field effects are negligible, in a rare example 
that is consistent with constraints on the tilt (see £4.4). Lower left: two-field effects are 
negligible, but many-field effects are large. 
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4.3.2 The incidence of multifield effects 



Equipped with a few examples of the possible phenomenology, we now proceed to a statistical 
analysis of the incidence of multifield effects in our ensemble. First, in Fig. [12] we show the 



cumulative probability that multifield contributions to the scalar power have a given size. The 
upper curve corresponds to £ e xact/one = Inexact /"Pone — 1|- We find that approximately 30% of 
realizations have £ eX act/one > 01, corresponding to a 1% correction to the power, and 10% of 
realizations have ^exact/one > 1? corresponding to a 100% correction to the power. The tail of 
the distribution is significant even at very large enhancements. 
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Figure 12: The relative probability P(£ > A) of different levels of corrections of the scalar 
power. The upper curve shows corrections due to multifield effects, corresponding to ^exact/one = 
l^exact/^one — 1|- The lower curve shows corrections due to many-field effects, corresponding to 
^exact/two = l^exact/'^two — 1|- Note that A = 1 corresponds to a 100% correction, and the tails 
extend to very large enhancements. 



Next, in Fig. [13] we show the relationship between the multifield correction to the power 
spectrum and the total turn A8 of the inflationary trajectory. An important fact visible in 
Fig. 13 is that for any value of the enhancement in power due to multifield effects, V exact /Vone, 



there is a minimum amount of total turn A9 necessary to effect such an enhancement. For 
total turning A8 exceeding a threshold value A6+ rs 1.4, the enhancement in power can be very 
large. For practical purposes, the vertical range of Fig. 13 has been restricted; see Fig. 12 for 
an alternative representation of the incidence of large enhancements. 
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Figure 13: The ratio of the exact power to the one-field power, versus the total turn of the 
trajectory during inflation, for all realizations. Notice the clear threshold of total turning 
required to achieve a given ratio of multifield to single-field power. 



4.3.3 The incidence of many-field effects 

We now turn to characterizing the incidence of many-field effects, characterized by ^exact/two = 
Inexact /"Ptwo — 1|- We remind the reader that the two-field truncation retains the instantaneous 
adiabatic and first entropic modes, so that ^exact/two > indicates that additional entropic modes 
contribute to the perturbations. For practical reasons, although we computed the exact scalar 
power for all 18731 realizations yielding at least 66 e-folds of inflation, we computed the two-field 
power only in the effectively multifield ensemble, for which ^exact/one = Inexact /V one — 1| > 0.01. 
The logic is that many-field effects will be significant only if multifield effects are significant, so 
that we can compute the frequency of many-field effects in the full ensemble without computing 
"Ptwo in the effectively single-field subset, 
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Many-field effects are only slightly less common than multifield effects in general: many- 
field corrections exceed 1% (i.e., ^exact/two > 0.01) in 18% of all realizations, and exceed 100% 
in 6% of realizations. Among effectively multifield models, many-field effects are commonplace: 



^exact/two > 0.01 in 62% of models with ^exact/one > 0.01. The lower curve in Fig. 12 indicates 
the cumulative probability that many-field contributions to the scalar power have a given size, 
as measured by £ exa ct/two- Comparing to the upper curve, which corresponds to ^exact/one = 
\P exac t /V one — 1|, we see that when multifield effects are large, the likelihood of many-field 



2(1 



In principle, one can envision a scenario in which V e xact ~ V ne <C 'Ptwo, so that two-field and many- 



field effects are both present but cancel out in the exact power (see £4.3.4 for a discussion of this point). 
Correspondingly, some scenarios with ^exact/one < 0.01 could involve significant many-field effects. We have 
verified that this occurrence is very infrequent in our ensemble, occurring in less than 1% of realizations. 
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effects grows dramatically. On the other hand, large many-field effects and large two-field 
effects are rarely present in the same model. The anticorrelation between two-field and many- 



field effects is clearly visible in Fig. 14 , which primarily consists of two distinct populations, one 



with exclusively two-field effects (horizontal branch) and another with exclusively many-field 
effects (vertical branch). It is worth stressing the peculiarity of this latter population: based on 
a two-field description, one would wrongly conclude that these models have negligible multifield 
effects, while in fact only higher entropic modes affect the curvature perturbation. 
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Figure 14: Many-field and two-field effects, for effectively multifield models (see £2.3.4). Two 
distinct populations are visible: the vertical branch has negligible two-field effects and sig- 
nificant many-field effects, while the horizontal branch has significant two-field effects and 
negligible many-field effects. Points below the red line have destructive multifield effects, i.e. 

Pexact < V one , Cf. mil 



4.3.4 Destructive multifield effects 



An interesting phenomenon visible in Fig. 13 and Fig. 14 is destructive interference, in which 
Pk+i < 'Pk for some k. Destructive two-field effects are not uncommon, but one can see from 
Fig. 14 that the majority of cases with "P two < V ori e have P e xact > "Pone, so that the net effect 
of all entropic modes is to add power. However, a handful of examples with P eX act < T^one are 



visible in Fig. 13 



In certain approximate treatments of the evolution of the perturbations in multifield infla- 
tion, the addition of entropic modes can only increase the power spectrum of the curvature 
perturbation 71, and the destructive multifield effects that we observe cannot arise. We will 
therefore pause to explain why destructive multifield effects are not forbidden, and what im- 
plicit assumptions lead to the conclusion that destructive multifield effects cannot arise. This 
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discussion is somewhat decoupled from the rest of the paper and can be omitted in a first 
reading. 

Let us consider a two-field inflationary model for the sake of simplicity. In the exact de- 
scription, a possible quantization scheme consists of writing 

Tir = Klal + (7^)*(a* k )t + ^ + (7^)*(a s _ k )t , (28) 
where the creation and and annihilation operators satisfy 

ft, (Sftl = - k') , (29) 



with I, J e (cr,s). The connection to the method of numerical evolution described in § |2.3.3 



is as follows. The exact coupled equations of motion are solved twice, for two different initial 
conditions. In a first run, the adiabatic fluctuation begins in the Bunch-Davies vacuum, while 
the entropic perturbation is initially set to zero. The corresponding solution for TZ leads to 
TZ k . In a second run, the entropic fluctuation begins in the Bunch-Davies vacuum, while the 
adiabatic perturbation is initially set to zero. The corresponding solution for 1Z leads to TZ k . 
The total power spectrum of 1Z is then given by 

(oi^r^Lr |o) = 5(k + k') (\K a k \ 2 + \n s k \ 2 ) . (30) 

Thus, in the exact description, there are two contributions to the curvature perturbation that 
add in quadrature. 

In the one-field description, on the other hand, one simply writes 

Ur = K k ai + (K k y(a'L^, (31) 

and the entropic fluctuations are set to zero for all time, which leads to 

(0|^ ne ^k" e |0) = S(k + k')\n k \ 2 . (32) 



If 1Z k in equation (28) and lZ k in equation (31) were identical, then we would conclude that 



entropic perturbations affect the two-point function of the curvature perturbation exclusively 



through the term \1Z k \ 2 in equation (30), and this contribution is manifestly nonnegative. 



In fact, however, 1Z k in equation (28) and lZ k in equation (31) are quite different. In the 
evolution that determines 7Z k , the entropic fluctuations are initially set to zero, but are not 
forced to be zero for all time: entropic fluctuations can be generated through coupling with 
the adiabatic mode, and can then backreact on the adiabatic mode itself. Thus, it can happen 
that | lZ k | 2 < \TZk\ 2 , and in turn it is logically possible — though not necessarily common — to 
have \n a k \ 2 + \TZ s k \ 2 < \TZ k \ 2 , i.e. V two < P one . 

Another perspective on this effect comes from considering the coupling between the adiabatic 
and entropic perturbations during Hubble crossing. A convenient assumption that is often 
utilized (but not always explicitly invoked) in the literature is that ?]i < 1 during the few e- 
folds of Hubble crossing, so that the adiabatic and entropic perturbations evolve independently 
during that time. In this approximation, adiabatic and entropic perturbations are uncorrelated 
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soon after Hubble crossing, say at a time t* at which spatial gradients can already be neglected, 
so that one can write (again taking a two-field model for simplicity) 

#k° = + (^)*(«- k J t + KtfL + (nr(a s _ k y , (33) 



where a%.+ and af^ satisfy the relation (29), and in particular are independent. We stress that 



and af^ are different from a£ and a£ in equation (28): the former are defined such that 
the power spectrum of TZ\ W ° at t* coincides with |7££(£*)| 2 . Just as before, Vf k {t > £*) is 
computed by solving the full system of equations, setting the entropic perturbations to zero 
at t — while TZ s k (t > t+) is computed by solving the full system of equations, setting the 
adiabatic perturbations to zero at t = t+. The power spectrum is then obtained as in equation 



(30). However, the assumption of decoupling during Hubble crossing results in one crucial 
difference, as we now explain. As is well known, although the entropic fluctuation can source 
the adiabatic fluctuation on super-Hubble scales, the adiabatic fluctuation does not source 
the entropic fluctuation in this regime [75]. Hence, if the entropic fluctuation and its time 
derivative are set to zero at time t+, they will remain zero, and no backreaction on TZ k (t > £*) 
is possible. Correspondingly, if multifield effects are neglected until a time t+ after which the 
adiabatic mode cannot source the entropic mode — which is precisely what occurs when the 
slow turn approximation is made during Hubble crossing — then TZ% coincides with the one-field 



quantity IZf. in (31). Because of the second term in (30), one then always obtains that, under 
these approximations, V two > V one . This corresponds, for example, to the prediction of the 
zeroth-order transfer functions formalism. 

As should be clear from the discussion above, a necessary condition for destructive multifield 
effects is bending of the trajectory around the time of Hubble crossing. As a consistency check, 
we have verified that in our realizations with destructive multifield effects, rj± always has a large 
peak at Hubble crossing. However, bending at Hubble exit is not a sufficient condition: a large 
number of our realizations have rj± > 1 during Hubble crossing but do not display destructive 
multifield effects. 

We should point out that although P exac t <^ V ori e occurs in our ensemble, V exact <C V na \ ve 
does not: when V exact <^ V one , we typically find V naive <^ V one and V exact = 0(V na \ ve ). Even so, 
examples with V exact < V na \ ve are common: 61% of all models have 0.63 < V exact /'P na \y e < 1. 
However, none of the realizations with V exact < V orie or V exact < V na \ ve has a scalar tilt in the 



observationally allowed region (see £4.4) 



Finally, we remark that the amplitude of primordial gravitational waves, as measured by r, 
can be slightly larger than what the naive estimate r = 16e would suggest. In several examples 
we indeed find VrfP exact > 16e, with Vj denoting the tensor power spectrum. It would be 
interesting to understand if it is therefore possible to violate the consistency relation r < —8n-f 
[371 176] . but to determine this requires computing nj- very precisely, which is challenging in 
models with r < 1. 

4.4 Constraints from scale invariance 

The results described in the preceding sections refer to the ensemble of realizations of inflation 
that give rise to 66 or more e-folds of expansion followed by a hybrid exit. Such a cosmic history 
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serves to solve the horizon, flatness, and monopole problems, but is not necessarily consistent 
with observations of the CMB temperature anisotropies. In this section we study multifield cor- 



rections to the tilt (£4.4.1), and then describe the restricted ensemble of realizations consistent 



with WMAP7 constraints on the tilt ( §4.4.21 ) . 



4.4.1 Multifield contributions to the tilt 

We begin by characterizing multifield contributions to the tilt n s of the scalar power spectrum. 



In Fig. 15 we show the correlation between the multifield correction to the tilt, defined as 
5n s = nl* act — n" aive , and the naive tilt n" aive . In realizations with a distinctly blue naive spectrum, 
inflation is occurring well above the inflection point, and transients from the approach to the 
inflection point have not necessarily decayed at the time that the CMB exits the Hubble radius 



(see { 2.2). Correspondingly, there is a larger likelihood of multifield effects, as evident in Fig. 15 
We observe that multifield effects typically shift the spectrum toward scale invariance, but the 
size of the effect is rarely large enough to produce a model in the observational window denoted 
by the orange band. There are a few very interesting cases with large multifield effects that 



have tilts 21 consistent with WMAP7 at the 2a level, which we will discuss in more detail in 



i 4.4.2 



The reader might object that Sn s = n^ xact — n" aive does not exclusively represent multifield 
effects: if n°" e ^ n" 3 ' ve then 5n s contains a contribution from strictly single-field slow-roll- 
violating effects. Slow-roll-violating effects, for effectively single-field trajectories, were already 
described in £3.3 (see Fig. [7]). We have verified that for the full ensemble, including realizations 



with large multifield effects, slow roll violations generally make small (5n^ R < 0.05) positive 



contributions to n s , as in the single-field case of { 3.3 In contrast, the multifield effects described 



above make a much larger negative contribution to n s . 

4.4.2 Realizations consistent with constraints on the tilt 

In this section we present results for the restricted ensemble of realizations that are consistent 
with WMAP7 constraints on the tilt of the scalar power spectrum, which constitute 21 % of our 
full ensemble. The overall amplitude of the scalar power spectrum can be adjusted by changing 



the scale _ of the potential, and we will assume that this has been done. We do not incorporate 
constraints on the running, but it could be interesting to do so. 

The consequence of imposing constraints on the tilt is very striking, and easily understood 



by recalling the single-field inflection point model described in §2.2[ A single-field inflection 
point model is compatible with constraints on the tilt only if the total number of e-folds is 
N e > 120, so that the observed CMB exits the Hubble radius when the inflaton is at or below 
the inflection point, at which point the spectrum becomes red. 

It follows that for any realization consistent with observations, either (i) N e > 120, and the 
tilt can be well-approximated by the single-field slow roll result, or (ii) the number of e-folds is 



21 We have not determined whether the running of the scalar power spectrum further constrains these models. 

22 In the underlying microphysical model, this adjustment corresponds to changing an inflaton-independent 
contribution to the vacuum energy. The magnitude of the vacuum energy could in principle be correlated with 
other terms in the potential, but modeling such a correlation is beyond the scope of the present work. 
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Figure 15: The difference between the exact and naive results for the tilt n s , versus the 



naive tilt, for effectively multifield realizations (see £2.3.4). Notice that multifield contributions 



typically shift the spectrum toward scale invariance, and the more blue the naive spectrum, 
the larger the scatter of 8n s . The orange band corresponds to the window allowed at 2a by 
WMAP7. 



constrained only by the solution of the horizon problem, and the tilt must differ substantially 
from the single-field slow roll result. Although case (ii) can occur, it is quite rare: multifield 
effects do tend to redden the spectrum, as noted in §4.4. 1[ but given how blue the naive spectrum 
is, the multifield effect is typically not large enough to bring the exact spectrum into agreement 
with observations (corresponding to points falling inside the orange band in Fig. 15). In case 
(i), which describes 98.7% of realizations consistent with observations, the 60 or more e-folds of 
inflation that precede the exit of the observed CMB generally suffice to damp out all transient 
effects from the approach to the inflection point (this is reflected in Fig. 16). Thus, although 
turning trajectories that imprint multifield effects in the curvature perturbations are present 
in a reasonable fraction (roughly 30%) of realizations of inflation, they only represent 1.3% of 
models with an observationally allowed tilt, in which these turns are generally complete long 
before the CMB exits the Hubble radius. 

For the same reason, violations of slow roll are rare in realizations consistent with observa- 
tions: 99% of models with an allowable tilt have —0.1 < m^/H 2 < —0.01. Similarly, trajectories 
with substantial bending are uncommon once constraints on the tilt are imposed. In the left 
panel of Fig. 16 we show r]± at Hubble crossing for the observationally allowed ensemble, which 
is to be contrasted with the left panel of Fig. 10 (In order to show the structure at small r]±, 
we have omitted a small number of points with rj± ~ 0.1.) The right panel of Fig. 16 shows 
the total turn of the trajectory, to be contrasted with the right panel of Fig. 10 Note that a 
handful of realizations allowed by constraints on the tilt do have significant bending. 

It is interesting to understand the characteristics of the realizations that have important 
multifield effects and are also consistent with constraints on the tilt. Perhaps surprisingly, these 
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realizations do not display extraordinary background behavior. Of course, all such models have 
non-negligible bending of the trajectory around or after Hubble crossing, but the maximum 
value of t]j_ along the inflationary evolution need not be greater than 0(0.1) to generate sig- 
nificant multifield effects. In fact, a majority display only a moderate degree of bending, 
?7_l ~ (9(0.05), albeit sustained over ten or more e-folds. One distinctive difference compared 
to the general population is a substantially smaller value of the first entropic mass, ml , at the 
time of Hubble crossing. Whereas m? Si ranges from to 10 in the full ensemble, typical values 
are of order 0.1 in observationally allowed realizations with important multifield effects. We 



show an example of such a realization in Fig. 17 
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Figure 16: Left panel: the slow roll parameter rj± evaluated at Hubble crossing, versus the total 
number of e-folds, N e , for models with tilt n s consistent with observations. Right panel: the 
total turn of the trajectory, from six e-folds before Hubble crossing until the end of inflation, 
A8 = j™£_ & dN e ri±, versus the total number of e-folds, N e , for models with tilt n s consistent 
with observations. 



4.5 Non-Gaussianities 

Although so far we have only discussed the scalar power spectrum, a major motivation for the 
study of multifield models is the prospect of a detectable primordial bispectrum. The D3-brane 
inflation scenario we have considered involves six fields that interact via Planck-suppressed 



couplings, leading to masses of order H for each field, ^ Our setting therefore corresponds 
to a microphysical realization of quasi-single-field inflation, introduced by Chen and Wang in 
[77] . Quasi-single-field inflation occupies the middle ground between single-field inflation, in 
which any scalar fields other than the inflaton have masses m 3> H, and the simplest models 
of multifield inflation, in which all the fluctuating fields have masses m <C H. 

A striking feature of quasi-single-field inflation is the possibility of large non-Gaussianity, 
given apparently modest cubic couplings in the Lagrangian, and modest rates of turning, as 



23 Recall from j |3.l| that this is only a parametric estimate, and there is considerable additional structure in 
the scalar mass spectrum. 
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Figure 17: An example in which multifield effects — specifically, two-field effects — are large, 
and the power spectrum is consistent with observations. This realization also appears in the 



lower right panel of Fig. 1 1 



measured by rj±. Specifically, there is a contribution to the bispectrum scaling as [77] 



1 V, 



NL 
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'IX 



(34) 



where V SlSlSl denotes the third derivative of the potential with respect to the instantaneous 
first entropic direction. Evidently, if V^ lSlSl is not very small compared to H then /jvz, can be 
large, even for perturbatively small T]±. 

We should point out that several important simplifying assumptions made in [77] are not 
necessarily applicable in our setting. In particular, rj± and the first entropic mass squared, m 2 , 
were assumed to be constant during inflation. Moreover, the authors of [77] considered a regime 
in which the multifield effects give subdominant contributions to the amplitude of the power 
spectrum (see also |78j). Finally, it was assumed in [77] that r/± was small enough to be used 
as a perturbative expansion parameter. This last point implies that the operator identified as 
leading in [77] may not give the dominant contribution to f^L in our context. In the following, 
we will obtain a rough picture of the prospects for non-Gaussianity in our ensemble using the 



estimate (34). Going beyond the approximations above is an interesting problem for the future. 



Although quasi-single-field inflation is a plausible proposal at the level of field theory, it 
is reasonable to ask whether the particular assumptions required to achieve detectable non- 
Gaussianity are natural in the effective theories arising from well-motivated microscopic con- 
structions. For clarity, we distinguish two properties of the effective theory: first, that typical 
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scalars — including the entropic fluctuations — have masses of order H; and second, that 



the cubic couplings V S1S1S1 are not too smal compared to H. There is very ample and 



long-established evidence for the first property in theories with spontaneously broken super- 
symmetry, as recently explained in detail in |3], and as noted above this is borne out in detail 
in our ensemble. However, the scale of the cubic couplings is more subtle, as we now explain. 

Consider a D3-brane in a conifold region of a de Sitter compactification of string theory. 
Using a spurion superfield X to represent the source of supersymmetry breaking, so that \F X \ 2 = 
3H 2 M 2 , the scalar potential for the D3-brane includes contributions of the form [3] 

V = J2ca / dto (|) ' = E 3c A ff2M P (|) ' /a(*) » (35) 

where A is an operator dimension, ca is a Wilson coefficient that is expected to be of order unity, 
and cf) = y/T^r, with T3 the D3-brane tension, is the canonically normalized field representing 
radial motion. Here /a(^) is a (known) function of the dimensionless angular coordinates 
0i> 0i, &2, 02, ip, collectively denoted 

It will be critical to relate the \& to the canonically normalized angular fields, which we 
denote by tp a , a = 1 ... 5. The metric on the angular manifold T 1,1 is not diagonal in the 
#1, 0i, #2, 02, "0 basis, but one can perform a coordinate transformation to obtain independent 
dimensionless fields, which we may write as \l/ a , a = 1 ... 5. The kinetic term for the ^f a takes 
the form C D J2 a=1 2 (<9^ a ) 2 , so at a given radial location <p = 0*, the canonically normalized 
angular fields are given by ip a = 0*^ a - 

The cubic coupling for three angular fields therefore takes the schematic form 

W~X> a g(A)"V(*), (36) 

where / A 6c denotes the corresponding derivative of /, which is generally of order unity. We 
observe that differentiating with respect to canonically normalized angular fields has introduced 
the factor (M p /0*) 3 . As many operators with A < 3 are present in the D3-brane Lagrangian, 
and the inflationary inflection point occurs at 0* <C M p , the cubic couplings of the angular fields 
are therefore parametrically large compared to ||-. This is worth emphasizing: if all fields had 
entered the Lagrangian on the same footing, and the scalar potential took the schematic for^ 1 ' 1 



V = H 2 M 2 P(^/M p , . . . <p N /M p ) , (37) 

with P a general polynomial of N fields <pi/M p , . . . , (p^/M p , with Taylor coefficients of order 
unity, then the cubic couplings would scale as V" ~ H 2 /M p . In the conifold (or any other 



24 For V SlSlS1 ^ H, there is a risk of large radiative corrections to the masses. We thank Daniel Baumann 
and Daniel Green for instructive discussions of this point. 

25 This estimate relies on the assumption that the fields are moduli whose masses and cubic couplings vanish 
before supersymmetry breaking. For D3-branes in the conifold, and more generally for moduli in string com- 
pactifications, this assumption is appropriate, but in more general settings larger cubic couplings need not be 
unnatural. We thank Daniel Baumann for helpful comments on this point. 
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Calabi-Yau cone) we find a potentially dramatic enhancement compared to this 'democratic' 
estimate. 

The discussion so far has only addressed the couplings of the angular fields. The cubic 
couplings of the radial field <fi include terms of the form 

H 2 / j. \ a-a 



^~c A A(A-l)(A-2)-^j . (38) 
Thus, if the A were all integers, we would have < H 2 /M p , which is too small to contribute 



to detectable non-Gaussianity (see (40) below). However, the operator dimensions are not all 
integers in the field theory dual to the conifold. Inflection point inflation in this model typically 
results from a cancellation between a term with A = 3/2 and the well-known conformal coupling 
term with A = 2 (see [H] for a discussion), which can occur at small values of 0. Differentiation 
of the term with A = 3/2 gives rise to a cubic coupling for (f> that is of the same order as the 
cubic couplings of the angular fields. 

To summarize, in a natural effective theory describing moduli that obtain their potential 
from gravitational-strength coupling to a source of supersymmetry breaking, in which the scalar 



potential involves a generic polynomial of the form (37), involving only integer powers of the 



fields, with all fields entering on the same footing, the cubic couplings scale as V" ~ H /M p , 
and are far too small to generate detectable non-Gaussianity. We have seen that D-brane 
inflation in the conifold necessarily modifies this simple picture, in two ways: the operator 
dimensions are not all integers, and the radial and angular fields play different roles in the 
potential. The result is that the cubic couplings of all six fields are parametrically of order 

H 2 / A X A-3 



V"~Y,TT TT • (39) 
Having addressed the parametric scalings of the cubic couplings, we turn to examining the 



range of numerical values that the quantities appearing in equation (34) attain in our ensemble 



As a measure of the scale H, we use e: equation (34) can then be rewritten as 



~ 10v^^%I ■ (40) 

The typical value is e ~ 10~ 12 , and there is a tail toward somewhat larger values: 7% of 
realizations have e > 10~ 8 , though we found no examples with e > 10~ 6 . The typical size of 
7]x is ~ 10~ 3 , while rj± < 1 in 99% of realizations, and we found no example with ?7_|_ > 3 26 
The typical scale of the cubic coupling of the first entropic mode in our ensemble is V S1S1S1 ~ 
10 3 if 2 /M p , and as a conservative estimate we take V SlSlSl < 10 4 H 2 /M p . Assembling these 
results, we find that /nl > 10 in less than 1% x 7% = .07% of all realizations. 

The results above apply to the full ensemble, but non-Gaussianity is far less likely in the 
ensemble consistent with constraints on the tilt: 97 % of observationally allowed models have 



26 These statements refer to the value of rj± averaged from one e-fold before Hubble exit until one e-fold after 
Hubble exit, which can be considerably smaller than the maximum value of r]± during this period. 
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< 10~ 3 , and 99% have e < 1CT 13 , so that there is a negligible prospect of detectable non- 
Gaussianity in such models. 

It is straightforward to extend our analysis to the quartic couplings, and correspondingly 
to estimate the overall amplitude t NL of the trispectrum. We rewrite the result [77] 

tsL ~ max (± ^V^^rA , (41) 

where V^ 1S1S1S1 denotes the fourth derivative of the potential with respect to the instantaneous 
first entropic direction, as 

tNL ~ 100 e r/_|_ max I I — I , — — M. (42) 

We find that typically V S1S1S1S1 < 10 8 H 2 /Mp, from which we deduce that there is a very limited 
possibility of a detectable trispectrum (t NL > 1000) in the full ensemble, but the corresponding 
probability is negligibly small in the ensemble consistent with constraints on the tilt. 

In conclusion, our results suggest that significant non-Gaussianity (/atl ^ 10 and/or tNL > 
1000) on scales far outside the present Hubble radius is possible, albeit somewhat rare, in 
warped D-brane inflation; but correspondingly large non-Gaussianity is very rare on observable 
angular scales. Even so, a dedicated study of non-Gaussianity in D-brane inflation — or in 
some other well-motivated microphysical realization of quasi-single-field inflation — would be 
worthwhile. 



5 Conclusions 

We have determined the scalar power spectrum that results when inflation takes place in a 
random six-field potential. The potential in question governs the motion of a D3-brane in a 
conifold region of a string compactification, and was derived and extensively studied in [5]. 
The signal property for this work is not the string theory provenance of the inflaton action, 
but merely the fact that the action is natural: the terms in the potential correspond to Planck- 
suppressed operators with coefficients of order unity, and the masses are of order H, as expected 
for moduli with gravitational-strength couplings to a source of supersymmetry breaking. For 
this reason, we believe our analysis is representative of a general class of multifield potentials 
that are natural in the Wilsonian sense, even though aspects of the conifold geometry do 
influence our results. 

The essence of our approach to computing the primordial perturbations was the comparison 
between numerical integration of the exact equations for the linearized perturbations, making no 
approximation, and an array of truncated models omitting one or more of the entropic modes. 
This allowed us to determine the extent to which various approximations and truncations - 
including the slow roll approximation, the single-field truncation, and the two-field truncation 
- capture the physics of a generic realization of inflation. As we made no slow roll or slow turn 
approximation, we were able to characterize phenomena that are common in our ensemble, but 
more rarely seen in analytic treatments. 
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We began by determining the spectrum of scalar masses using a random matrix model, 
building on [68], and demonstrated excellent agreement with simulations of the full potential. 
We found that at the time of Hubble exit, one entropic mode is typically light enough to 
fluctuate, but before the end of inflation an adiabatic limit is reached, and one can predict the 
late-time curvature perturbation without modeling the details of reheating. 

Our results for the perturbations are usefully divided into characterizations of the ensemble 
of realizations of inflation that give rise to at least 60 e-folds of inflation, but are otherwise 
unconstrained, and characterizations of the ensemble of realizations of inflation that give rise 
to at least 60 e-folds of inflation, and are also consistent with observations of the CMB temper- 
ature anisotropies: conditioning on the approximate scale-invariance required by observations 
drastically changes the outcome. 

For the ensemble of inflationary models that solve the horizon problem but are not required 
to give nearly scale-invariant density perturbations, multifield effects are often significant: mul- 
tifield corrections to the spectrum are at least at the 1% level in roughly 30% of realizations, 
and exceed 100% in 10% of realizations. We also found that many-field contributions to the 
perturbations — by which we mean effects that cannot be described by the two-field trunca- 
tion, which retains only the instantaneous adiabatic and first entropic modes — are similarly 
common: many-field corrections exceed 1% in 18% of realizations, and exceed 100% in 6% of 
realizations. Most models with substantial multifield effects have either large two-field effects 
or large many-field effects, but not both. Finally, we observed that the exact scalar power can 
be smaller than the power calculated in the single-field truncation: this is a consequence of 
trajectories that turn quickly at the time of Hubble crossing. 

Turning now to the restricted ensemble of realizations that are consistent with the WMAP7 
constraints on the tilt of the scalar power spectrum, corresponding to 21% of the total ensemble, 
we find a qualitatively different picture. When inflation occurs near an approximate inflection 
point, as it does in the class of potentials studied here, the tilt depends on whether the observed 
CMB exits the Hubble radius before or after the inflaton descends past the inflection point. 
When the 60 e-fold mark occurs above the inflection point, the curvature of the potential is 
positive and the spectrum computed in the single-field slow roll approximation is blue, which 
is inconsistent with observations. Therefore, the allowable potentials are those in which 60 or 
more e-folds occur below the inflection point. In such potentials, there is a prolonged phase 
of inflation before the observed CMB exits the Hubble radius: namely, the inflation occurring 
above the inflection point. This phase generally suffices to damp out all transients from the 
onset of inflation, including the curving trajectories that can give rise to multifield effects in the 
perturbations. For this simple reason, although multifield contributions to the perturbations 
are commonplace in random inflection point models, for more than 98% of models with an 
observationally allowed tilt these contributions arise on angular scales that are far outside 
the present Hubble radius. We conclude that for inflection point models of the sort studied 
here, multifield effects in the observable perturbations are uncommon in realizations consistent 
with present limits on scale-invariance. We stress that the multifield effects themselves are 
not responsible for the problematic deviations from scale invariance: in fact, the spectrum 
calculated in the single-field approximation is blue, and multifield effects shift it toward scale 
invariance. Instead, multifield effects generally arise from transients from the onset of inflation, 
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and at the onset of inflation at an inflection point, the dominant single- field component of the 
spectrum is often, but not always, unacceptably blue. 

The ensemble we have studied provides a concrete microphysical realization of quasi-single- 
field inflation. The couplings in the inflaton action arise from Planck-suppressed operators, and 
the cubic and quartic interactions of the entropic modes, although enhanced by contributions 
from an operator with dimension A = 3/2, as well as by the intrinsic asymmetry between 
the radial and angular directions of the conifold, are only occasionally sufficiently large to 
produce detectable non-Gaussianity along the lines envisioned in [77] . More generally, the 
suppression of transients resulting from constraints on the tilt ensures, as explained above, that 
in a majority of realizations consistent with observations of the spectrum, the single-field slow 
roll approximation is valid. Of course, the perturbations are Gaussian to good approximation in 
all such cases. Understanding the generality of this finding in well-motivated effective theories 
is an important problem for the future. 

We expect that our methods, and some aspects of our findings, have broad applicability. Our 
present understanding of string compactifications suggests considering ensembles of effective 
theories with numerous moduli and spontaneously broken supersymmetry. In this work we 
have characterized the cosmological signatures of one such ensemble. We have seen that when 
inflation occurs at an inflection point in a many-field potential, then even though it is common 
for multiple fields to fluctuate, and ultimately to contribute to the curvature perturbations, 
these contributions are most often on unobservably large angular scales. The signatures of the 
model are therefore often indistinguishable from those of a single-field inflection point scenario. 
Nevertheless, it remains to be seen which other universality classes of multifield inflation may 
exist, and it is reasonable to anticipate very different conclusions when the number of fields is 
large. 
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